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We calculated the von Neumann entanglement entropy and the Schmidt number of one dimen- 
tional (ID) cluster states and showed that these are useful measures to estimate entanglement caused 
by delocalization of clusters. We analyze system size dependence of these entanglement measures 
in the linear-chain na states given by Tohsaki-Horiuchi-Schuck-Ropke wave functions for ID cluster 
gas states. We show that the Schmidt number is an almost equivalent measures to the von Neumann 
entanglement entropy when the delocalization of clusters occurs in the entire system but it shows 
different behaviors in a partially delocalized state containing localized clusters and delocalized ones. 

It means that the Renyi-2 entanglement entropy, which relates to the Schmidt number, is found 
to be almost equivalent to the von Neumann entanglement entropy for the full delocalized cluster 
system but it is less sensitive to the partially delocalized cluster system than the von Neumann en¬ 
tanglement entropy. We also propose a new entanglement measure which has a generalized form of 
the Schmidt number. Sensitivity of these measures of entanglement to the delocalization of clusters 
in low-density regions was discussed. 


I. INTRODUCTION 

Nuclear many-body systems are self-bound systems of four species of Fermions, spin up and down protons and 
neutrons. There, a variety of phenomena arise originating in many-body correlations. One of the remarkable spatial 
correlations is ’’cluster” which is a subunit composed of strongly correlating nucleons. A typical cluster in nuclear 
systems is an a cluster, which is a composite particle consisting of four nucleons. If there is no correlation between 
nucleons in a nucleus, all nucleons behave as independent particles and the nucleus is an uncorrelated state with a 
clear Fermi surface written by a single Slater determinant wave function. However, in reality, correlations between 
nucleons are rather strong and a clusters are often formed at the surface in particular in light nuclei. Once a clusters 
are formed, delocalization of clusters occurs to gain the kinetic energy of center of mass motion (c.m.m.) of clusters 
in some situations. The delocalization of clusters involves many-body correlations beyond a Slater determinant. 

In realistic nuclear systems, a degree of the (de)localization of clusters depends on the competition (balance) of the 
kinetic energy and potential energy of clusters and is strongly affected also by the Pauli blocking between nucleons 
in clusters and a core nucleus. The delocalization limit of a clusters is an a cluster gas state, where all a clusters 
move almost freely like a gas. Such a cluster gas has been predicted to appear in the second 0 + state of 12 C 0 , 0 - To 
describe cluster states of delocalized a clusters, a new type of cluster wave function, the so-called “Tohsaki-Horiuchi- 
Schuck-Ropke” (THSR) wave function, has been introduced by Tohsaki et al. [2]. The THSR wave function is 
essentially based on a clusters in a common Gaussian orbit having a range of the system size, and suitable to describe 
general cluster gas states of n a clusters. Indeed, it has been shown that 8 Be(0^) and 12 C(0^) can be described well 
by the THSR wave functions of 2 a and 3a, respectively. Since the optimized THSR wave functions for these states 
have much larger ranges of the Gaussian orbit of clusters than the cluster size, 8 Be(0+) and 12 C(C >2 ) are interpreted 
as gas-like cluster states of 2a and 3a HI- 

The THSR wave function has been extended to apply to 20 Ne, and found to be able to describe also 16 0+a states 
in 20 Ne @,0. Recently, Suhara et al. have proposed that this concept of the a-cluster gas is applicable also to one 
dimension (ID) cluster motion in linear-chain na structures [§]. They have proposed the 1D-THSR wave functions 
and shown that the 1D-THSR wave functions with the optimized Gaussian ranges can describe well the exact solutions 
of the linear-chain 3a and 4a states. Their result somewhat contradicts to the the conventional picture of the linear- 
chain na structures that spatially localized a clusters are arranged in ID with certain intervals |9|. Even though their 
model restricted in ID is not enough to settle the problem of stability of the linear-chain states in realistic nuclear 
systems in 3D, their work provides a new picture of ID cluster states and may lead to a understanding of cluster 
phenomena in nuclear many-body systems. Moreover, the ID cluster state is an academically interesting problem of 
quantum many-Fermion systems. 

To distinguish between localization and delocalization of composite particles (clusters) in microscopic wave functions 
of Fermion (nucleon) systems, one should carefully consider the antisymmetrization effect of nucleons between clusters. 
When clusters largely overlap with each other, motion of clusters is strongly affected by the Pauli blocking between 
nucleons in other clusters. As a result of Pauli blocking effect, when the system size is as small as or smaller than the 
cluster size, clusters can not move freely and the system becomes equivalent to a localized cluster state. In the case 
of a large system size, where the overlap between clusters is small, clusters can move rather freely. It means that the 
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delocalization occurs not in high-density regions but in low-density regions. Indeed, we have investigated the cluster 
motion of an a cluster around the 16 0 and shown that the delocalization of the cluster occurs in a long tail part far 
from the core 0 - Also in the ID na states, Suhara et al. have shown peak structures in the density distributions of 
the linear-chan 3a and 4a states suggesting partial localization at high density inner regions even in the ID gases of 
3a and 4a. 

To estimate many-body correlations caused by the delocalization of clusters we need a measure of correlations that is 
free from the antisymmetrization effect. There are several entanglement measures that are essentially based on density 
matrix and can be used to estimate many-body correlations in quantum many-body systems. The Schmidt number 
11 is one of the measures, and the von Neumann entanglement entropy jl2| is another measure. The former has been 
proposed by Grobe et al. in 1994 EH, and applied to measure entanglement, i.e., many-body correlations in such 
systems as atomic and nuclear systems (see, for example, Refs. jl3l - fl5j and references therein). The latter has been 
introduced by Bennett et al. in 1996 [12|, and widely used in various fields such as condensed matter and quantum 
held theory (see Refs. [l6i - [20j and references therein). In the previous paper, we have calculated the entanglement 
entropy for the one-body density matrix in ID cluster states to measure the entanglement (correlations) caused by the 
delocalization of clusters^and found that the delocalization occurs in low-density regions but is relatively suppressed 
in high density regions [21]. 

Let us consider uncorrelated and correlated states of Fermion systems. An uncorrelated state is given by a Slater 
determinant wave function, for which the one-body density matrix p is a projector p 2 = p in the single-particle 
Hilbert space. It means that the eigen values pi of one-body density matrix satisfies p 2 = pi, i.e., pi = 1 or 0. pi is the 
occupation probability of the single-particle basis that diagonalizes the one-body density matrix and it is nothing but 
the Schmidt coefficients in the Schmidt decomposition for the one-body density matrix. Both of the Schmidt number 
and entanglement entropy are entanglement measures which can estimate how the Schmidt coefficients deviates from 
the condition p 2 = pi for uncorrelated states. However, these measures have different dependences on the one-body 
density matrix, i.e., the eigen values pi. 

In this paper, we calculate the entanglement entropy and the Schmidt number of ID cluster states and make it 
clear whether these entanglement measures are useful to investigate the delocalization of clusters. We analyze the 
system size dependence of these measures in ID a-cluster states given by the 1D-THSR wave functions. We show 
that the entanglement entropy and the Schmidt number are almost equivalent measures when the delocalization of 
clusters occurs in the entire system but they show different behaviors in partially delocalized cluster states composed 
of localized clusters and delocalized ones. We also propose a new entanglement measure which has a generalized form 
of the Schmidt number. 

This paper is organized as follows. We describe the entanglement measures defined by the one-body density matrix 
in Section[TTl Section[HT]discusses properties of these measures of ideal states in a toy model. In section HVl we calculate 
these measures in ID nuclear systems of a clusters and discuss sensitivity of these measures to the delocalization of 
clusters. The paper concludes with a summary in section [Vj 


II. MEASURES OF ENTANGLEMENT 

The entanglement entropy and Schmidt number are the measures of entanglement in quantum many-body systems, 
which are defined by the density matrix. In this section, we describe these measures and also propose a new measure 
by extending the Schmidt number. We also define spatial distributions of these entanglement measures. 


A. One-body density matrix 

In the present paper, we use only the one-body density matrix, which has been often discussed in nuclear systems 
p3 |. though more general density matrices are used to define measures of entanglement. For a wave function |’F^) 
of an A-particle state, the matrix element of the one-body density in an arbitrary basis is given as 

p pq = (^ [A) \c\ c p\¥ A) ), (1) 

where cj and c p are the creation and annihilation operators, respectively. The one-body density matrix is regarded 
as the matrix element of the one-body density operator p^ for the wave function db" 4 ), 

fa\ =^\p)ppq(q\- 

pq 


( 2 ) 
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The one-body density matrix can be diagonalized by a unitary transformation of single-particle basis 

(D^pD)ui = pi5w, (3) 

a l = Di'icji , (4) 

v 

where 

pi = (’F (A) |a]'az|'I' ( ' 4) }, (5) 

0 < pz < 1 (6) 

is the eigen value of the one-body density matrix and means the occupation probability of the single-particle state 
l in the wave function T ( ‘ 4 b pi corresponds to so-called Schmidt coefficients in the Schmidt decomposition for the 
one-body density matrix. The trace of the one-body density matrix p equals to the particle number A as 

A = Tr P = 5>- (7) 

i 

Note that the normalization of the one-body density matrix is not a unit but is the particle number A. In the 
present paper, the entanglement entropy and the Schmidt number are defined by thus defined one-body density 
matrix normalized as Trp = A. 

If a wave function |5'^ 4 ^) is a Slater determinant, p; = 1 for occupied single-particle states and pi = 0 for unoccupied 
single-particle states. It means that the one-body density operator p^ satisfies p\ = p^ and is a projector in the 
single-particle Hilbert space if is a non-entangled state given by a Slater determinant wave function. 


B. Measures of entanglement 

1. Entanglement entropy 

The von Neumann entanglement entropy has been introduced bv Bennett et al. and proved to be an entanglement 
measure in quantum many-body systems (see, for instance, Refs. [15l. fl8j and references therein). The entanglement 
entropy is defined by the von Neumann entropy of one of the reduced density matrices and called “von Neumann 
entanglement entropy”, which we call the “entanglement entropy” unless otherwise noted. In the present paper, 
we consider the entanglement entropy only for the one-body density matrix of Fermion systems. The entanglement 
entropy that is defined by the one-body density matrix is given as, 

S = -Trp log p = pi log p L . (8) 

l 

The entanglement entropy is zero if a wave function |’I'(' 4 )} is a Slater determinant, because pi = 1 for occupied 
single-particle states and pi = 0 for unoccupied single-particle states. It means that a system has non-zero positive 
value of the entanglement entropy only if the system contains many-body correlations beyond a Slater determinant, 
i.e., if the system is entangled. 


2. Schmidt number 


Another entanglement measure is the so-called Schmidt number, which has been introduced by Grobe et al. and 
used to measure many-body correlations in atomic and nuclear physics fl3Hl5l |. The Schmidt number K is defined as, 


K = 


A 

eTpF 


A 

Trp2’ 


( 9 ) 


which estimates the number of states involved in the Schmidt decomposition. I\ equals one if a wave function IH/^) 
is a Slater determinant, and K is greater than one for entangled states. In analogy to the entanglement entropy, 
which is generated by the entanglement, it is useful to consider the quantity K — 1, 

^ i Ei(pi-pj) 

Ei pi Ei pi 


K - 1 = 


( 10 ) 
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It is clear that K — 1 = 0 only for non-entangled states because pf = pi, i.e., p\ = py is satisfied only if a state is 
given by a Slater determinant. It means that a non-zero positive value of K — 1 is generated by entanglement, that is, 
many-body correlations beyond a Slater determinant. One of the merits of the K number is that it is given by Trp 2 
which can be calculated by the matrix element p pq in an arbitrary basis without the diagonalization of the one-body 
density matrix. 

We should comment that the logarithm of the Schmidt number is nothing but the Renyi entanglement entropy of 
order 2 (Renyi-2 entanglement entropy) for the one-body density matrix, 

log K = S 2 Renyi , (11) 

Sf enyi = ( 12 ) 

It is known that, in the £ —> 1 limit, the Renyi entanglement entropy becomes equal to the von Neumann entanglement 
entropy. In this paper, we discuss the Schmidt number instead of Renyi-2 entropy though they are equivalent 
entanglement measures. 


3. Extension of Schmidt number 


In the present paper, we propose an entanglement measure which is regarded as a generalized version of the Schmidt 
number. As shown in Eq. m, the origin of non-zero contributions in K — 1 is partially occupied single-particle state 
with 0 < pi < 1. Ignoring the total scaling factor 1 /EzP?j the contribution W{pi) (the weight function) of a single¬ 
particle state with the occupation probability 0 < p < 1 in K — 1 has the pi dependence, W(pi) = pi — p\, which is 
different from the weight function W(pi) = —pi log pi in the entanglement entropy S. The former has relatively small 
weights for single-particle states with small occupation probability pi < 1/2 compared with the latter although the 
occupied pi = 1 and unoccupied pi = 0 single-particle states have no weight in both cases (see the p dependence of 
the weight functions in Fig. [[]). We here introduce an alternative weight function, W(pi) = p/( 1 — pi), and define an 
entanglement measure K 1 by extending the entanglement measure K as follows, 


IF, 

K 7 - 1 


EiPl Trp 7 

p| 1+7) Trp( 1+7 ) ’ 

E Pi) 

E/p| 1+7) ’ 


(13) 

(14) 


where the parameter 7 is a positive constant. In the case of 7 = 1, the Af 7 number becomes consistent with the 
Schmidt number K. Similarly to the known entanglement measures (S and K — 1), AT 7 — 1 equals to zero only for 
a Slater determinant, whereas a non-zero positive value of Kj — 1 is generated by entanglement, that is, many-body 
correlations beyond a Slater determinant. Note that the diagonalization of the density matrix is required to obtain 
the Ky number for a non-integer 7 differently from the K number. In this paper, we choose 7 = log 2 which gives 
the weight function W(pi) = p° g2 {l — pi ) having pi dependence similar to W(pi) = —pi log pi for the entanglement 
entropy in 0.001 < p < 1 (see Fig. [TJ). 

It should be commented that log K 7 is given by the Renyi entanglement entropy of order 7 and 1 + 7 as, 


log K 1 = (l- 7 )S Renyi + 7 E R “ yi - 


(15) 


C. Spatial distributions of entanglement measures 

To investigate the spatial distribution of the ’’important single-particle states” that contribute to the non-zero 
entanglement entropy, we have defined the spatial distribution s(r) of the entanglement entropy, which have been 
introduced in the previous paper [EJ as, 

S = ^(-pilogpi) = f s(r )dr, (16) 

1 

s(r) = ^2 (-pi log pi)(f*( r )M r )- 

l 


(17) 
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FIG. 1: Weight functions W(p) = —plogp, p(l — p), and p log2 (l — p) for S, K, and K\ og 2 , respectively. 

Here the factor —pjlogp; is contribution of the single-particle state |Z) in S, and 0 z *(r)</>/(r) means the density 
distribution in |Z) and it is normalized as f 4>i(r)(j)i(r)dr = 1. Therefore, the distribution s(r) reflects spatial 
distributions of the important single-particle states |Z) that contribute to the total entanglement entropy, whereas it 
is hardly affected by almost occupied single-particle states having pi « 1. The expression of S with s(r) is analogous 
to that of the particle number with the density distribution, 

A = [ ^( r ) dr > (18) 

l J 

p ( r ) = (19) 

i 

Note that the distribution s(r) is not quantity determined only by local information at the position r. It is different 
from the density distribution p(r) which is determined only by the local information. 

We also define the distributions n(jr) and k 7 (r) for K — 1 and K 1 — 1, respectively, to see spatial distributions of 
the ” important single-particle states” in total amount of the measures, 

EM-pj) 

E iPf 

= J n(r)dr , ( 20 ) 

= H X-' r 2 ( Pl ~ Pf)<t>l(. r )M r ) 

I El Pi 

- ~ p'lWiMMr), ( 21 ) 

l 

EjtfXj - pi ) 

Ei p ] +1 

J K 7 (r)dr, ( 22 ) 

X! W \+ 7 p7(! - Pi)<l>i(r)Mr)- 
i Ev Pp 


K -1 


n(r) 


and 

AT, - 1 = 


K 7 (r) = 


(23) 
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Similarly to s(r), these distributions, n(r) and Ky(r), show spatial distributions of the important single-particle states 
|Z) that contribute to the total amount of K — 1 and Ky — 1, respectively. Note that these distributions are not local 
quantities. 


III. S, K, AND K-, OF IDEAL STATES IN A TOY MODEL 

In this section, we discuss behaviors of S , K, and Ah, for correlated (entangled) states composed of delocalized 
clusters in a toy model. Here, 7 in the Ky number is assumed to be 0 < 7 < 1. 

Let us consider n/ species of particles. For instance, n/ = 2 for spin up and down Fermions, and n/ = 4 for spin up 
and down protons and neutrons. We use the label a for species of Fermions such as a = p f, P 4-, n f, n 4- for nuclear 
systems. We consider an A-body system containing the same number n = A/rif of a particles. We assume that the 
system is symmetric for the exchange of species, that is, single-particle orbitals are occupied by all species of particles 
with an equal weight. Then the density matrix is diagonal with respect to < 7 , and all species of particles have the same 
occupation probability p a i = pi independent from a. We consider the one-body density matrix and operator in the 
reduced space with the dimension n and define the entanglement measures, S, K, and Ky , for a species of particles 
by using the reduced matrix of the one-body density. The total entanglement entropy Sfotai and the total numbers 
A’totai and A 7)to tai are given by the measures for each species as Sfotai = n/S, AT to tai = K, and A' 7jto tai = AT. In this 
paper, we discuss S, AT, and Ky for a species of particles. 

For simplicity, A particles are assumed to stay on sites in a space. The number of available sites (single-particle 
states) is m and we use the label kj for the jth site. Let us first consider a system of A = rif particles. If an A-body 
state is an ideal state of independent particles, the wave function can be written by a simple product of single-particle 
wave functions 


^(1,2, • • •,«/) = (24) 

0 (*) = c(k)<f> k (i). (25) 

k—ki 

For instance, c(k) is constant as c(k) = 1/y/rn for a free gas state. For such an non-entangled state, S = 0, AT — 1 = 0, 
and Ky — 1 = 0 because the one-body density operator is given as p ^ = |V’)('0| and satisfies p\ = p Another 
example is a “localized cluster” system of a cluster, where all particles are localized at one site kj to form a composite 
particle (a cluster) at kj. The wave function is given as 

n f 

*(1,2.n/) = n^( 0 - (26) 

2=1 

This wave function for a localized clusters also has zero measures, S = 0, AT — 1 = 0, and Ky — 1 = 0, because the 
one-body density operator is given as p^ = | <j>kj){<i>kj I and satisfies again p\ = p^. Namely, the localized cluster wave 
function is a non-entangled state. 

We consider a state of a delocalized cluster in a strong correlation limit, 

{ n f n f n f ^ 

n ^ (*)+n 0** (*)+••• n 0*- w [ > ( 2? ) 

2=1 2=1 2=1 J 

where the composite particle composed of A = n/ particles moves freely in the entire space with an equal probability 
■A. This is a highly entangled (strongly correlated) state, where, if a particle is observed at a certain site, all other 
particles are always observed at the same site. This is a strong coupling limit and an example of a delocalized cluster. 
The one-body density operator, 

771 1 

fa =^2~\ k o)(kj\, (28) 

i=i m 

corresponds to the Schmidt decomposition with the common Schmidt coefficients, 1 /to. We get S = log to, K = to, 
and Ky = to. In general, if eigen values of pi are constant pi = 1 /my (1 = 1,... ,mv) for a given number (my) of 
states and pi = 0 for l > my, we get 


S 

I< 

Ky 


log my, 


TOy, 

TOy. 


e S = TOy 


(29) 

(30) 

(31) 
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It indicates that e s , K , and Ab, equal to the number my of states involved in the Schmidt decomposition, my is 
regarded as the effective volume size. Note that, for a Slater determinant, e s , AT, and A ' 7 equal to 1 indicating that 
the effective volume size is one which can not be decomposed. 

Next we consider an ^4-particle system of n = A/rif clusters. Here n is the number of clusters (composite particles) 
formed by rif constituent Fermions. For a state of n localized clusters at n sites k = kj 1 ,, kj n , the wave function 
is given as, 


{ rt—l 

II <t>k H ( hn f + 1) • • • (t> kjl ( hn f + n f ) 
h =0 

= {n\)~ nf/2 A \^ kjh (1) • • • <t> kH (n f ) . <j> kjn (A-n f + 1) • • • (j) kjn (A)|. (32) 

Since the occupation probability is pi = 1 for occupied single-particle states kj x ,..., kj n and pi = 0 for unoccupied 
single-particle states, the system has the zero value of the measures, S = 0, K — 1 = 0, and K 1 — 1 = 0. 

Let us consider a state of n clusters in the delocalized limit where all clusters are delocalized and move freely in 
a volume size my like a gas. For this state of delocalized clusters, the occupation probability is pi = n/my for 
l = 1,... ,my. my should not be less than n because of the Pauli principle of n Fermions. The measures of this 
delocalized cluster system are, 


S 

I\ 

AT 


m v 

n log-, 

n 

my 

5 

n 

my 

n 


e S/n = 

n 


(33) 

(34) 

(35) 


It indicates that K and A ' 7 are consistent with e s ^ n . ne s ^ n , nK, and nK 1 equal to the number my of the states 
involved in the Schmidt decomposition and they estimate the effective volume size of the delocalization of clusters. 
For the case of my = n = A/nj, all my single-particle states are completely occupied by A particles and clusters 
can not move at all. The state is equivalent to the localized cluster, and it has zero value, S = 0, K — 1 = 0, and 
K 1 — 1 = 0, of entanglement measures. In case my is larger than n, the measures, S , K — 1, and K 7 — 1, become 
positive indicating that the delocalization of clusters occurs and the system becomes an entangled state. In the case 
that my is much larger than n, the system corresponds to a low-density cluster gas and it is a highly entangled state. 

Finally, we consider the case of a partial delocalization that n — 1 clusters are localized to form a core and only 
the last cluster is delocalized. This situation corresponds to an a cluster moving almost freely around a core nucleus. 
In this partially localized case of a delocalized cluster around a core, the occupation probability is pi = n/my for 
1 = 1,..., my and pi = 1 for the n — 1 single-particle states occupied by constituent particles of the core. Then, the 
measures of this partially delocalized system are 


S = log my, 

K _ nmy 


» s / n = m\' 


AT 7 = 


(n - 1 )my + 1 ’ 


(n — l)r 


my 


(■n — 1 )my + 1 


In the low-density limit of the large my, we get 


K 

AT 7 


n 

n — 1 ’ 


m. 
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(36) 

(37) 

(38) 


(39) 

(40) 


Let us compare the fully delocalized cluster state and the partially delocalized cluster state in the large my limit 
(the large volume size limit). The former state corresponds to a dilute cluster gas, where all clusters are delocalized 
moving freely like a gas, and is a highly entangled state. For the fully delocalized cluster state, all the entanglement 
measures are sensitive to my as given in Eqs. HMD, m , and (1351) . That is, e s ^ n , AT, and K 1 equal to my In and 
equivalently good indicators to measure the delocalization of clusters. However, for the latter case of the partial 
delocalization, which corresponds to a delocalized cluster around a core, e s / ra , K, and A ' 7 are not equivalent but 
show different dependences on my. As clearly shown in Eqs. (ESI). ESI), and (TO . e s /n anc j ^ increases as the 
my increases, but K becomes constant and does not depend on my in the large my limit. It indicates that S and 






Kj can be useful measures sensitive to the partial delocalization in subsystem, but K is insensitive to the partial 
delocalization. The reason for the insensitivity of K is the significant contribution from fully occupied single-particle 
states with pi = 1 in the denominator of the definition of the K number, which makes the contribution from the 
delocalized part minor. From another point of view, it is found that K can be a good probe to clarify whether the 
delocalization of clusters occurs in the entire system or not. my dependences of e s ^ n and if 7 are the powers of 1 /n 
and 1 — 7 , respectively. More generally, for the partially delocalized system that is composed of n g delocalized clusters 
and no localized clusters, the my dependences of e s ^ n and K 1 in the large my limit are the powers oin g /n and 1 — 7 , 
respectively. Here, n g and no are the numbers of delocalized and localized clusters, respectively, and no + n g = n. 
This means that, in the case n g /n > 1 — 7 of a small fraction of delocalized clusters, K-y is more sensitive to the 
delocalization of clusters than e' 5 /™, whereas, in the case of n g /n « 1 — 7 , the K 7 number has the my dependence 
similar to e s ^ n . 


IV. APPLICATION TO ID NUCLEAR SYSTEMS OF a CLUSTERS 

In the present paper, we use the delocalized cluster wave functions in ID for the linear-chain na states which are 
investigated in the previous paper [2l[. We also adopt the a + (2a) wave functions for a state of an a cluster around 
a 2 a core. We analyze the entanglement measures , S , K, and ifi og 2 , and also their spatial distributions, and discuss 
the system size dependence of these measures. 


A. Localized and delocalized a cluster wave functions in ID 


We here briefly explain the adopted model wave functions for (de)localized cluster states in ID. More details of the 
model wave functions are described in the previous paper [21 .] . We consider intrinsic wave functions of the linear-chain 
na-cluster states aligned to the x axis. It means that the (de)localization of a clusters are defined for a-cluster motion 
along the x axis. 

For a localized na-cluster wave function, we use the BB wave function [23] as, 


(Rl, . . . , R n ) 

- [V& ■ ■ ■ V’flJ > 

(41) 

r Ri 

= 4 >( R i Xpt ( l >< R i Xpi4 ,( R i Xn-f(l> ( R i Xni', 

(42) 

(/) os 

= (nb 2 )- 3/4 exp — ^2 (? Ri) 2 ■ 

(43) 


'i/’fl. is the four-nucleon wave function of the ith a cluster expressed by the (Os) 4 harmonic oscillator (ho) shell-model 
configuration localized around the spatial position R,. % is the spin-isospin part and </> 0s is the spatial part of the 
single-particle wave function. For ID cluster states, the position parameter R, is set to be i?, = (7?.,;, 0, 0), and the ID 
BB wave function is expressed as $§b(I?i, ... , R n ). The parameter b for the a-cluster size is chosen to be b = 1.376 
fin same as in Ref. [s]. Note that a single BB wave function is a localized cluster wave function written by a Slater 
determinant of single-particle wave functions and it has exactly zero values of the measures, 5 = 0, K — 1 = 0, and 
K. y — 1 = 0. General wave functions for ID na systems can be written by linear combination of BB wave functions 

‘I’mV R\ ./Gi. 

For a delocalized cluster wave function, we use the 1D-THSR wave functions of na , which have been proposed by 
Suhara et al. to describe the linear-chain 3a- and 4a-cluster states in 12 C and 16 0 systems jsj. The 1D-THSR wave 
functions are given by linear combination of BB wave functions with a Gaussian weight as, 



If the antisymmetrization is ignored, the 4>™Q_ THgR (/3) expresses the na state where all a clusters are confined in the y 
and z directions while they move in the x direction in the Gaussian orbit with the range parameter f3. (3 corresponds 
to the system size of the ID na state. In case of the system size /3 is as small as or smaller than the a-cluster size 
&, the 1D-THSR wave function is approximately equivalent to a localized cluster wave function given by a Slater 
determinant because of the antisymmetrization effect. As (3 increases, the delocalization of a clusters occurs. When 
/? is large enough compared with the a-cluster size b, the system goes to a dilute ID a-cluster gas where n a clusters 
move almost freely like a gas in the x direction. 
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In the practical calculation, the Ri integration is approximated by summation on mesh points in a finite box as done 
in the previous paper. For 2a, 3 a, and 4 a systems, we make a correction of the total c.m.m. to eliminate a possible 
artifact from /3 dependence in the total c.m.m. as described in the previous paper. It means that the 1D-THSR wave 
function of la without the c.m.m. correction expresses a system of an a cluster bound in an external field, which is 
not a realistic state of an isolate nucleus, whereas those of 2a, 3a, and 4a with the c.m.m. correction correspond to 
self bound na systems with linear-chain structures predicted in nuclear states such as excited states of 12 C and 16 0. 
In this paper, we discuss the entanglement measures for la with no c.m.m. correction (nc) and those for 2a, 3a, and 
4a with the c.m.m. correlation. We also show the enganlement entropy for 2a with no c.m.m. correction, just for 
comparison. 

We also use the 1D-THSR wave function of a + (2a) for the case of an a cluster around the 2a core, 

$ id- 2 thsr(^) = y'dRiexp|-^-|$| < |(I?i,I ?2 = +£,R 3 = ~e), (45) 

where the 2a core is located at the origin and an a cluster is distributed around the core with a Gaussian weight. 
When f} is large enough compared with the a-cluster size b, the wave function describes a delocalized a cluster around 
the 2a core at the origin. This wave function is associated with the partially delocalized cluster wave function in the 
toy model discussed previously. We use a small value of e = 0.02 fm to describe the 2a core almost equivalent to the 
h.o. (0s) 4 (0pa ;) 4 configuration. This wave function has been originally introduced in previous paper to describe the 
a+ 16 0 cluster states in 20 Ne. 

We calculate numerically the one-body density matrices for these ID cluster wave functions, and obtain the mea¬ 
sures, S, I\, and /\i og 2 from the eigen values pi of the one-body density matrices. The detailed method of the practical 
calculation is described in the previous paper. Spatial distributions, s(r), k (r), and K 7 (r) and the density distribution 
p(r) are integrated out along y and z axes, and we get distributions, s(x), n(x), n-y(x), and p(x), projected onto the 
x axis. 


B. na cluster states 

We analyze S, K, and I \\ og 2 of the ID na states written by the 1D-THSR wave functions, 4 >id_ thsr (/ 3). Figure 
[ 2 ] shows the system size dependence of the entanglement entropy of the la state with no c.m.m. correction (la nc ), 
and 2a, 3a, and 4a states with the c.m.m. correction, as well as that of a 2a state with no c.m.m. correction (2a nc ), 
S/n, e s / n , and n(e s ^ n — 1) are plotted as functions of the dimensionless system size (3/b. In the /3/b = 0 limit, the 
entanglement entropy equals zero, because the 1D-THSR wave functions are equivalent to localized na cluster wave 
functions in this limit. As the (3/b increases, the delocalization of clusters occurs and the entanglement entropy is 
generated. e s ^ n for the la nc state increases almost linearly to (3/b as expected from a naive picture that a cluster 
moves freely in a finite volume, which is decomposed into my ~ (3/b states by the quantum decoherence in the one- 
body density matrix. Similarly to the la nc state, S for 2a, 3a, and 4a states increases with the increase of the system 
size, indicating that the entanglement entropy is generated as the delocalization of clusters is enhanced. e s ^ n for 2 a, 
3a, and 4a states also shows almost linear dependence to (3/b. The n(e s ^ n — 1) plots in Figs. [IJc) and (21(d) show a 
good correspondence of the (3 dependence of the entanglement entropy between la nc and 2a nc , and that between 2a, 
3a, and 4a states. The values of n(e s ^ n — 1) for 2a, 3a, and 4a are relatively smaller compared with those for la nc 
and 2a nc because of the c.m.m correction performed for the 2a, 3a, and 4a states. 

Let us compare other measures, K and A'i og 2 , with the entanglement entropy. Figure |3] shows the system size 
dependence of e s ^ n , K, and ATi og 2 of la nc , 2a, 3a, and 4a states. K and A'i og 2 shows the /3/b dependence quite 
similar to that of e s ^ n except for global normalization factors. This result indicates that the entanglement entropy, 
K, and K\ og 2 , can be approximately equivalent measures for the cluster delocalization of na states. It is consistent 
with the naive expectation from the analysis of the delocalized cluster states in the toy model discussed in the 
previous section. It means that ne s ^ n , nK, and nK \ og 2 estimate the number of the states involved in the Schmidt 
decomposition as shown in Eqs. (1331) . (1341) . and (1351) . 

We also compare the spatial distributions of these measures, s(a;), n(x), and n\ og 2 (x) of la nc in Fig. [IJ The 
density distribution is also shown. As described previously, s(x), n(x), and n\ og 2 (x) reflect spatial distributions of the 
important single-particle states that give non-zero contributions to total measures S, K — 1, Ai og 2 — 1, respectively. 
Note that the shape, in particular, the spatial broadness of distributions is of importance, but their global scaling 
(normalization) is not so meaningful. It is found that s(cc) and K\ og 2 (x) show quite similar distributions to each other. 
They are more broadly distributed than the density distribution indicating that S and A'i og 2 — 1 are generated in low- 
density regions but relatively suppressed in high-density regions. Differently from s(x ) and Ki og 2 (a : ), the enhancement 
in low-density regions and the suppression in high-density regions are relatively weak in n(x). As a result, the spatial 
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extent of k(x) is not as remarkable as that of s(x ) and Ki og 2 (^)- This difference in the spatial distributions of S , K — 1, 
R) og 2 — 1 comes from the different weight functions W (p) , namely, S and K\ og 2 — 1 are relatively sensitive to single¬ 
particle states with low occupation probability compared with K — 1 as already shown in Fig. [l] The distributions 
s(x), k(x), and ni OS 2 (x) for the 2a, 3a, and 4a states are shown in Fig. [5] as well as density distribution. Similarly to 
the la„ c state, s(x) and Ki og 2 (x) are more broadly distributed than the density distribution and also slightly broader 
than k(x), indicating again that S and K \ og 2 — 1 are generated in low-density regions but suppressed in high-density 
regions. 

In the present result, we find that the S, K and K i og 2 can be useful measures to estimate the entanglement caused by 
the delocalization of clusters in the ID na states given by the 1D-THSR wave functions. As the system size increases, 
the delocalization of clusters develops and non-zero S , K — 1, and ATi og 2 — 1 are generated. In the spatial distributions 
of these entanglement measures, significant contributions come from low-density regions than high-density regions. 
Quantitatively, the Schmidt number ( K ) is less sensitive to low-occupation probability single-particle states than 
the entanglement entropy ( S ) and the Ad og 2 number resulting in the less broad distribution of k(x) than s(x) and 
K\ 0 g 2 (x), which are remarkably broader than the density distribution. 



FIG. 2: System size /3/b dependence of the entanglement entropy S and e s ^ n in the 1D-THSR wave functions of na. n(e s ^ n — 1) 
is also shown. The 1D-THSR wave functions with the total c.m.m. correction are used for 2a, 3a, and 4a systems, and that 
with no c.m.m. correction is used for the la system (la n c). The results of the 1D-THSR wave function with no c.m.m. 
correction for 2a system (2a nc ) are also shown for comparison. 


C. an a cluster around a 2a core 

We investigate S, K , and i\i og 2 of a+(2a) cluster wave functions. We use the 1D-THSR wave function 4>“ ] i ) (2 rrHSR(/3) 
in Eq. (H51) for the a + (2a) states associated with the partially delocalized cluster wave function. We also adopt the 
parity projected BB wave function with the fixed a-(2a) distance d as used in the previous paper as, 

Kf ah+ (d) = (l + P r )$l c h(Ri=d,R 2 = +e,R 3 = -e) 

= ^bb(-^i = R2 = T £ > R3 = ~ £ ) 

+$!bCRi = -d, R 2 = +e, Rs = -e), (46) 

where P r is the parity transformation operator. Note that $gg 2a ^’ + (d) is given by the linear combination of two 
Slater determinants. This corresponds to the symmetry restored state where the parity symmetry is broken in the 
intrinsic state because of the cluster development. 

Figure [5] shows S, K , and A'i og 2 as well as e S// " for the ID THSR wave function and the parity-projected BB wave 
function of the a + (2a) system. Looking at the result of the parity-projected BB wave function, we find that, as d 
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FIG. 3: Comparison of system size /3/b dependences between e 5 ^ n , K , and K \ og 2 in the 1D-THSR wave functions of not. 
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FIG. 4: Spatial distributions, s(ir), n(x), and Ki 0 g 2 (*), of S, K — 1, and iS'iog 2 — 1 in the 1D-THSR wave functions of la n c 
with P = 1 fm, 3 fm, and 5 fm. The corresponding dimensionless system sizes are (3/b = 0.73, 2.18, and 3.63. The density 
distribution p(x) is also shown. 


increases and the parity symmetry is broken in the intrinsic wave function, non-zero values of S', K — 1, and ivi og 2 are 
generated in the projected state even though the intrinsic wave function before the parity projection is the localized 
cluster wave function. When d is enough large, P r <&|b becomes independent to $bb> and we get S —»• log 2 = 0.693, 
K — 1 — > 1/5 = 0.2, K\ og2 — 1 — > 0.236 from the eigen values of the one-body density matrix, p 1 = p 2 = 1 and 
P 3 = Pi = 1/2. It means that non-zero values of S, K — 1, and Ri og 2 — 1 are generated by the symmetry breaking and 
restoration. Let us turn to the result of the 1D-THSR wave function. In the (3/b = 0 limit, S, K — 1, and K\ og2 — 1 
are zero. In the (3/b < 1 region, S, K, and Ki og2 increases rapidly with the increase of (3/b because of the symmetry 
breaking and restoration as seen in the parity-projected BB wave function. In the (3/b > 1, S, K , and Ad og 2 increases 
gradually as the system size increases indicating that the delocalization of cluster develops in this region. 

As discussed in the previous section, the Schmidt number, K , should be less sensitive to the delocalization in the 
partially delocalized cluster states. To see the sensitivity of e s ^ n , Ji, and Ki og2 to the delocalization, we show in 
Fig. Hi) the scaled measures, (e s / n — l)/0.26, ( K — l)/0.2, AsTi og 2 /0.236, in the 1D-THSR wave function of a + (2a), 
which are normalized to the values of the large d limit of the parity-projected BB wave function. As expected from the 
analysis of the simple toy model, the result in Fig. Hi) shows that K is not so sensitive to the delocalization of cluster 
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FIG. 5: Spatial distribution, s(x ), of the entanglement entropy S in the 1D-THSR wave functions of 2a, 3 a, and 4a with j3 = 2 
fm, 5 fm, and 8 fm. The corresponding dimensionless system sizes are fl/b = 1.45, 3.63, and 5.81. The density distribution 
p(x) is also shown. 
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FIG. 6 : Same as Fig. [5]but spatial distribution, «(#), for K — 1. 


in the a + (2a) state, whereas e s / n and K\ og 2 more strongly depend on the the system size in the /3/b > 1 region 
than K. The j3/b dependence of K \ og 2 is quite similar to that of e s ^ n , maybe, because of the accidental coincidence 

of the m|/" dependences, e s ^ n oc mj/" and AT 7 oc TOy -7 ^ f° r the partially delocalized cluster state discussed in the 
simple toy model as l/n = 1/3 and 1 — 7 = f — log 2 = 0.31. 

Figure |9] shows the distributions s(x), k(x), and Ki og 2 (^) of S, K — 1, and K 1 — 1 as well as the density distribution 
in the 1D-THSR wave functions of a + (2a) for /? = 2 fm and 5 fm. It is clear that distributions s(x), k(x), and 
Kiog 2 (a 0 are strongly suppressed in the |x| < 2 fm because of the Pauli blocking effect from the 2a core. In the /3 = 5 
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FIG. 7: Same as Fig. [3] but spatial distribution, K\ og 2 (x), for i^iog 2 — 1. 


fm case, the delocalization occurs, and S and K 7 — 1 is generated in the low-density regions in particular at the long 
tail part. Similarly to the distributions in na states, the distribution k(x) for K — 1 is not so enhanced in low-density 
regions as s(x) and Ki og 2 (z) because the Schmidt number K is relatively less sensitive to single-particle states with 
low occupation probability compared with the entanglement entropy and the K \ og 2 number. 
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FIG. 8 : (a)(b)(c)(d) System size /3/b dependence of S, e s ^ n , K, and A'i og 2 in the 1D-THSR wave function of a + (2a), and 
(e)(f)(g)(h) d/b dependence of S, e s ^ n , K, and Aj og 2 in the parity-projected BB wave function of a + (2a). (i) Scaled measures, 
( e s / n _ 1)/0.26, (K — l)/0.2, Ajog 2/0.236, in the 1D-THSR wave function of a + (2a), normalized to the large d limit values 
of the parity-projected BB wave function. 
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FIG. 9: Spatial distributions, s(x), k(x), and Ki og 2 (*), of S, K , and A'i og 2 in the 1D-THSR wave functions of a + (2a) with 
/3 = 2 fm and 5 fm. The corresponding dimensionless system sizes are f)/b = 1.45 and 3.63. The density distribution p(x) is 
also shown. 


V. SUMMARY 

We calculated the entanglement entropy ( S ) and the Schmidt number ( K ) defined by the one-body density matrix 
in the ID a-cluster states to measure the entanglement caused by the delocalization of clusters in nuclear systems. 
We also propose a new entanglement measure with a generalized form of the Schmidt number. 

For the delocalized cluster states of na given by the 1D-THSR wave functions, e s /", K , and Ad og 2 show good 
correspondence indicating that the entanglement entropy, the Schmidt number, and the K\ og 2 number are almost 
equivalent entanglement measures to estimate the delocalization of a ID gas of delocalized n a clusters. On the 
other hand, for the partially delocalized cluster state which contains a delocalized cluster and a core composed of 
localized clusters, the Schmidt number is not sensitive to the delocalization of the cluster around the core, whereas 
the entanglement entropy and the ATi og 2 number can be good indicators to measure the partial delocalization. In 
other words, the Schmidt number can be a good probe to clarify whether the delocalization occurs for all clusters 
in the entire system or not, owing to the insensitivity to the partial delocalization. We should comment that the 

. . nRenyi 

Schmidt number corresponds to the the Renyi-2 entanglement entropy; K = e * 2 . It means that the equivalence 

in the fully delocalized cluster states and difference in the partially delocalized cluster states between e s ^ n and K are 
nothing but those between the von Neumann entanglement entropy and Renyi-2 entanglement entropy. 

In the present analysis of the ID a-cluster states, the /\i og 2 number shows similar features to the entanglement 
entropy. It indicates that the K\ og 2 number can be an alternative measure to the entanglement entropy to estimate 
the delocalization in both cases of the partially and fully delocalized cluster states. When the delocalized part is minor 
compared with the core part, the Ai og 2 number is more sensitive to the partial delocalization than the entanglement 
entropy, and hence it is a promising measure. 

We should point out that the calculation of the Schmidt number may be practically easier than those of the 
entanglement entropy and K 1 with a non-integer 7 because the Schmidt number is given by Trp 2 for which the 
diagonalization of the one-body density matrix is not required. This could be an advantage of the Schmidt number in 
large dimensional calculations because numerical errors might become a more serious problem in practical calculations 
of the entanglement entropy and the AT 7 number. 

Recently, comparison between von Neumann and Renyi-2 entanglement entropies have been discussed in various 
fields such as field theories j24j. The equivalence and difference between von Neumann and Renyi-2 entanglement 
entropies shown in the present study are general features and may be found in various quantum systems. The present 
study of entanglement measures in cluster wave functions of nuclear systems may shed light on study of entanglement 
(correlation) in general quantum systems. 
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